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Abstract 

Processes of heavy quark photoproduction at HERA energies are considered using 
the semihard {k± factorization) QCD approach with emphasis of the BFKL dynam- 
ics of gluon distributions. We investigate the dependences of the total cross section 
of heavy quark photoproduction and also pT, and rapidity distributions on different 
forms of the unintegrated gluon distribution. We present a comparision of the theo- 
retical results with available HI and ZEUS experimental data for charm and beauty 
photoproductions . 
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1 Introduction 



Recently, HI and ZEUS Collaborations have reported []T], |^ experimental data on the total 
ctoss section of inelastic open beauty photoproduction. A comparision of these results with 
NLO pQCD calculations shows that ones underestimate the cross section at HERA energies. 
Therefore, it would be certainly reasonable to try a different way. 

In the present note, we focus on the so called semihard approach |] (SHA), which we 
had applied earlier to open heavy quark ^ and J /'^ 01 photoproduction. 

At the HERA energies and beyond, the interaction dynamics is governed by the properties 
of parton distributions in the small x region. This domain is characterized by the double 
inequality s 3> /i^ — s ^ A^, which shows that the typical parton interaction scale /i (mass 
rric or px of heavy quark) is much higher than the QCD parameter A, but is much lower 
than the total c.m.s. energy ^/s. The situation is therefore classified as "semihard". 

The resummation i, i |, i of the terms [ln(/iVA2) a,]", [ln(/iVA2) ln(l/x) a,]" and 
[ln(l/x) cts]" in SHA results in the unintegrated parton distributions /x), which 

determine the probability to find a parton of type i carrying the longitudinal momentum 
fraction x and transverse momentum qx at the probing scale /i^. They obey the BFKL 



equation and reduce to the conventional parton densities Fj(a;, /x') once the qx dependence 
is integrated out: 

/ ^i{x,ql,^?, Ql) dql = xFi{x,^l'', Ql). (1) 
Jo 

To calculate the cross section of a physical process, the unintegrated functions $j have 
to be convoluted with off-mass shell matrix elements corresponding to the relevant partonic 
subprocesses. In the off-mass shell matrix element the virtual gluon polarization tensor is 
taken in the form of the SHA prescription |P, 

=^ = p'p'x'/lqTl' = q^TQTnqT? (2) 

In ref. |^ was used phenominological parametrization for the unintegrated gluon distri- 
bution including an arbitrary normalization constant {K— factor), which was obtained from 
a fit to 66— pair production at the Tevatron 0. In our recent paper [0 we investigate the 
sensitivity of heavy quarkonium photoproduction to different gluon distributions. Special 
attention was given to the unintegrated gluon distributions obtained from BFKL evolution 
equation. In this paper we study the sensitivity of the total cross section of inelastic charm 
abd beauty photoproduction to these unintegrated gluon distributions. The outline of our 
paper is the following: In sect. 2, we give the formulas for the cross sections of heavy 
quark photoproduction in the SHA of QCD. Then, in sect. 3, we describe the unintegrated 
distributions which we use for our calculations. In sect. 4, we present the results of our 
calculations. Finally, in sect. 5, we give some conclusions. 



2 SHA QCD cross section for heavy quark photopro- 
duction 

We calculate the total and differential cross sections (the p±_ and rapidity distributions) 
of charm and beauty quark photoproduction via the photon-gluon fusion QCD subprocess 
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(Fig.l) in the framework of the SHA. 

First of all we take into account the transverse momentum of gluon (f2±, its the virtuality 
g| = — gfj^ and the alignment of its polarization vectors along its transverse momentum such 
as = g2±/i/ I ^2± I ii S El- define Sudakov variables of the process ep QQX 

(Fig.l): 

pi = aiPi + (3iP2 + piL P2 = a2-Pi + P2P2 + P2± 

qi = XiPi + qi± q2 = X2P2 + q2± (3) 

where 

22;\^2 22 22 

Pi=P2 = M , q^= gi^, q2 = q2±, 

Pi and p2 are 4-momenta of the heavy quarks, qi is 4-momentum of the photon, q2 is 4- 
momentum of the gluon, p2±; g2± are transverse 4-momenta of these ones. In the 

center of mass frame of colliding particles we can write Pi = {E, 0, 0, E), P2 = {E, 0, 0, —E), 
where E = P^ = P2 = and {P1P2) = s/2. Sudakov variables are expressed as 

follows: 

Mi_L , M2± , 

tti = ^^exp(|/J a2 = ^^exp(?/2) 

Pi = ^^exp(-?/i) /32 = ^^exp(-?/2), (4) 

where Mf2± = M'^ +Pi2±5 2/12 rapidities of heavy quarks, M is heavy quark mass. 
From conservation laws we can easly obtain the following conditions: 

qi± + q2± = Pi± + P2±, xi = ai + a2, X2 = Pi + P2 (5) 
The differential cross section of heavy quark photoproduction has the following form 



(Ppi± J TT 167r2(sa;2)^a2 

The matrix element M for a subprocess 7(7* qq depends on the virtuality of the gluon 
and differs from the one of the usual parton model. For the square of this matrix element 
we used the following form M: 



where s, t, u are usual Mandelstam variables of partonic subprocess 7(7* qq. 



3 Unintegrated gluon distribution 

In this paper we used the different parametrizations for the unintegrated gluon distribution. 
First, as in the publication 0, we used the following phenomenological parametrization 
(LRSS-parametrization) Q: 

qr ') = $ 077#^(1 " ^) Qt (8) 
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where 



with q^{x) = qi + exp(3.56y'/n(xo/x)), = 2GeV^, A = 56MeV, Xq = 1/3. The value 
of the parameter q^ix) can be considered as a new typical transverse momentum of the 
partons in the parton cascade which leads to natural infrared cut-off in semihard approach. 
The normalization factor $o of the structure function $(x,gT^) was obtained in 0, where 
beauty production at CDF energy was described, $o = 0.97 mb. 

The effective gluon distribution xG{x,fi'^), which was obtained from eq. (7)-(8), and eq. 
(1) increases at not very low x(0.01 < a; < 0.15) as 

xG(x,/i2) ~ x"^, (10) 

where A ^ 0.5 corresponds to the QCD pomeron singularity given by summation of leading 
logarithmic contributions {agin-)"' This increase continues up to a; = xq, where xq is a 
solution of the equation q^ixo) = yU^. In the region x < xq, there is saturation of the gluon 
distribution function: xG{x,fi'^) ~ $o/^^- 

A second parametrization is based on the numerical solution of the BFKL evolution 
equations |1T[ (RS-parametrization). The solution has the following form [|11|]: 



f)= [a2 + a3(%) + (%)^ 

«2 + 0-3 + a4 q^ q^ 



e + ln(l/x) + X 

a^x'"' In (gVa8)](l + aiix)(l - a;)"9+"ioi'^(9'/''B), (11) 

where 

C =[^' if < qo{x), ,^2) 

" \qo{x)/q^, iiq^>qo{x). 



All parameters (see |Tl|]) (ai — an, a, (3 and e) were found by minimization of the differ- 
ences between left hand and right-hand of the BFKL-type equation for unintegrated gluon 
distribution $(x,g^), Ql =4 GeV^. 

Finally we also use the results of a BFKL-like parameterization of the unintegrated 
gluon distribution $(x, /i^), according to the prescription given in |T^. The proposed 
method lies upon a straightforward perturbative solution of the BFKL equation where the 
coUinear gluon density xG(x,/i^) from the standard GRV set |12] is used as the boundary 



condition in the integral form (1). Technically, the unintegrated gluon density is calculated 
as a convolution of collinear gluon density G(x, /i^) with universal weight factors ||12|| : 



^x, ql, fx') = f gir^, ql, /i^) - G(-, /i^) dr^, (13) 



where 



a. 



g{v,qT,f^n = ^ M'2sJasHyv)Hf^yq'T)), ql < /^', (i4) 
vqr 



Of course LRSS and RS parametrizations are BFKL - type too. 
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(y I 

where Jo and Jq stand for Bessel functions (of real and imaginary arguments, respectively), 
and as = 3as/7f. The parameter as is connected with the Pomeron trajectory intercept: 
A = (5^4 In 2 in the LO and A = Q;s41n2 — Na'^ in the NLO approximations, respectively, 
where iV ~ 18 The latter value of A have dramatic consequences for high energy 



phenomenology. In particular it leads to negative values for physical cross sections |115 
However some resummation procedures proposed in the last years lead to positive value of 
A(~ 0.2 — 0.3) [l^. Therefore in our calculations with (13) we used only the solution of 



LO BFKL equation and considered A as free parameter varying it from 0.166 to 0.53. 

The presence of the two different parameters, /i^ and q'^, in eq.(12) for unintegrated gluon 
distribution $(a;, gf^, /i^) refers to the fact that the evolution of parton densities proceeds in 
two steps. First the DGLAP scheme [|1^ is applied to evolve the structure function from 
to within the collinear approximation. After that eqs. (12)-(14) are used to develop the 
parton transverse momenta q^ in correspondence with BFKL evolution [|l2 . 



This approach was used for the description of pr spectrum of D* meson electroproduction 
at HERA 1 19 1 and J/\l/ photoproduction [0, where in the first case for Pomeron intercept 
parameter was obtained the value A = 0.35. However for the total cross section of inelastic 
J/\E' photoproduction A = 0.53 is more preferable. 



4 Results of calculations 

The calculations of the heavy quark photoproduction cross section in the SHA have been 
made accoding to eqns. (6) and (7) for ^ > Qo GeV^ and for qx^ < q^ GeV^ we set 
Ig^l = in the matrix element of process and use one of the standard partom model (SPM). 
The choice of the critical value of parameter qx"^ = q^ = 1 — 2 GeV^ is determined by the 
requirement of the small value of as(gT^) in the region ^ > 1 — 2 GeV^, where in fact 
as(gT') <O.26 0. 

The results of our calculations for the cross sections of cc and hh photoproduction pro- 
cesses are shown in Figs. 2 - 10. Fig. 2 shows the total cross section of the inelastic cc 
photoproduction at HERA and fixed target energies as a function of ^/s. The curves 3, 4 
and 5 correspond to the SHA calculations with the RS, LRSS and BFKL unintegrated gluon 
distributions at the Pomeron intercept A = 0.35, = 1.5 GeV and at values of = 4, 2 
and 1 GeV^ accordingly. One can see that the SHA curves 4 and 5 describe the data better 
than the SPM ones (1 and 2). The RS parameterization of the unintegrated gluon distri- 
bution gives the total cross section more close to SPM result (obtained with GRV gluon 
density) because of very large value of = 4 GeV^. 

The LRSS parameterization of the unintegrated gluon distribution gives the cross section 
which stops to increase at energy y/s > 150 GeV because of saturation effects accounted for 
by the parameterizations (8) and (9). We would like to stress that account for saturation 

^The substitution = 1 GeV^ instead 2 GeV'^ essentially influences on the dependence of cross sections 
from the BFKL parameterization of unintegrated gluon distribution. 
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effects in the region x < 10 ^ 20], does not contradict the HI and ZEUS experimental 
data |2T|,g. 

The BFKL parameterization (curve 5) describes the HI and ZEUS data very well at 
rric = 1.3 GeV (Fig. 3) The curves 3, 4 and 5 obtained in the SHA show more rapid growth 
with energy than those obtained in the SPM with GRV [1^ (curve 2) and MT (curve 1) 
parameterizations. 

In Fig. 4 we show the dependence of the total cross section from A parameter for the BFKL 
parametrization (13) - (15) at rric = 1.3 GeV. The dependence is very weak in contrast with 
the dependence of the total inelastic J/^ photoproduction cross section 0. 

The results of our calculations for the px and rapidity y {y = yl) distributions of the c 
quarks at \/s = 200 GeV are shown in Figs. 5 and 6. The SHA curves 4 and 5 in Fig. 5 for 
the pt distribution are higher than those of the SPM. It is effect of the SHA broadening of 
Pt spectra due to extra transverse momentum of the initial gluons 0. The SHA rapidity 
distributions (Fig. 6) have more similar behaviour in wide rapidity region. 

Fig. 7 shows the total cross section of inelastic bb photoproduction at HERA energies 
as a function of ^/s. The curves 3, 4 and 5 correspond to the SHA calculations with the 
RS, LRSS and BFKL unintegrated gluon distributions at the Pomeron intercept A = 0.35, 
TUb = 4.75 GeV and at values of = 4, 2 and 1 GeV^ accordingly. The HI [|ll and EMC |2§ 
experimental data are described by the LRSS parameterization (curve 4) and the BFKL 
parameterization (curve 5 in Fig. 8) but at small value of rrif, = 4.25 GeV only for the last 
parameterization. 

The results of the SHA calculations for the pt and rapidity distributions of the b quarks 
at ^/s = 200 GeV are shown in Figs. 9 and 10. The differences between SHA curves and 
SPM ones are more pronounced than in the c-quark case. 



5 Conclusions 

We considered the process of inelastic heavy quark photoproduction at HERA in the frame- 
work of the semihard QCD approach with emphasis on the BFKL dynamics of gluon distri- 
butions. We investigated total cross section of inelastic c— and 6— quark photoproduction 
and also pt and rapidity spectra as a function of different unintegrated gluon distributions. 
It is shown that the description of the c— quark inelastic cross section at HERA energies is 
achieved in the cases of the LRSS and BFKL parameterizations at A = 0.35 and at the 
values of rric = 1.5 and 1.3 GeV and at gg = 2 and 1 GeV^ accordingly. 

The cross section of inelastic 6— quark photoproduction at HERA is described by the 
LRSS parameterization at the value of = 4.75 GeV and A = 0.35 very well. The BFKL 
paremeterization with A = 0.35 describes the HI and EMC experimental data at ruf, = 4.25 
GeV only. 

The dependence of the total cross section of inelastic open c— and 6— quark photopro- 
duction on A parameter for the BFKL parameterization (13) - (15) is very weak in contrast 
with the dependence of the total inelastic J /"if photoproduction cross section. 

^ Close values for the parameter A were obtained rather in very different papers (see, for example, |25j 
-[p8[) and in the L3 experiment [ p9| . 
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Figure captions 

Fig. 1 QCD diagram for the open heavy quark electroproduction. 

Fig. 2 The total cross section of inelastic cc photoproduction as a function y/s 
at rric — 1.5 GeV. 

Fig. 3 The total cross section of inelastic cc photoproduction as a function ^/s 
at nic — 1.3 GeV. 

Fig. 4 The total cross section of inelastic cc photoproduction as a function ^/s 
at mc = 1.3 GeV and different values of A: 1 - A = 0.166, 2 - A = 0.35, 
3 - A = 0.53. 

Fig. 5 The distribution of inelastic cc photoproduction at ^/s —200 GeV. 

Fig. 6 The rapidity distribution of inelastic cc photoproduction at ^/s =200 GeV. 

Fig. 7 The total cross section of inelastic bb photoproduction as a function ^/s 
at nib = 4.75 GeV. 

Fig. 8 The total cross section of inelastic bb photoproduction as a function ^/s 
at rub = 4.25 GeV. 

Fig. 9 The p'^ distribution of inelastic bb photoproduction at ^/s —200 GeV. 
Fig. 10 The rapidity distribution of inelastic bb photoproduction ^/s =200 GeV. 
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Fig. 5. 
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Fig. 10. 
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